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A  conical  surface  roughness  model  applicable  to  particle  simulations  has  been  developed.  The 
model  has  been  experimentally  validated  for  channel  flows  using  helium  and  nitrogen  gases  at 
Reynolds  numbers  from  0.01  to  10  based  on  inlet  conditions.  To  efficiently  simulate  gas-surface 
interaction,  molecular  collisions  with  the  actual  rough  surface  are  simulated  by  collisions  with  a 
randomly  positioned  conical  hole  having  a  fixed  opening  angle.  This  model  requires  only  one 
surface  parameter,  average  surface  roughness  angle.  This  model  has  also  been  linked  to  the 
Cercignani-Lampis  scattering  kernel  as  a  required  reference  for  use  in  deterministic  kinetic  solvers. 
Experiments  were  conducted  on  transitional  flows  through  a  150  ptm  tall,  1  cm  wide,  1.5  cm  long 
microchannel  where  the  mean  free  path  is  on  the  order  of  the  roughness  size.  The  channel  walls 
were  made  of  silicon  with:  (i)  polished  smooth  surfaces,  (ii)  regular  triangular  roughness,  and  (iii) 
regular  square  roughness  with  characteristic  roughness  scales  of  <1  /im,  11  yim,  and  29  yim, 
respectively.  For  the  triangular  roughness,  mass  flow  reductions  ranged  from  6%  at  the  higher 
stagnation  pressures  tested  to  25%  at  the  lower  stagnation  pressures  tested  when  compared  to  the 
smooth  channel.  ©  2007  American  Institute  of  Physics.  [DOI:  10.1063/1.2775977] 


I.  INTRODUCTION 

Surface  structure  and  roughness  plays  a  role,  to  some 
extent,  in  all  flow-surface  interactions  regardless  of  speed  or 
scale.  Traditionally,  surface  roughness  is  considered  impor¬ 
tant  in  low  speed  microflows,  where  the  surface  to  volume 
ratio  is  large  and  boundary  conditions  dominate  the  flow.1 
Since  the  first  report  was  published,  almost  a  century  ago,"  a 
substantial  amount  of  research  on  the  impact  of  surface 
roughness  on  gas  flows  in  ducts  has  been  accomplished.  For 

instance,  the  mass  flow  degradation  due  to  surface  roughness 

3 

was  shown  numerically  and  experimentally  by  Davis  et  al: 
More  recently,  significant  effect  of  surface  roughness  on 
pressure  and  heat  transfer  in  circular  tubes  and  channels  was 
demonstrated  experimentally  in  Mala  and  Li4  and  Kandikar 
et  al?  In  the  experiments  by  Turner  et  al.6  however,  the 
surface  roughness  was  found  to  have  a  small  effect  on  gas 
pressures  inside  a  channel  for  several  pressure  ratios.  In  ad¬ 
dition  to  low  speed  flows,  surface  roughness  has  been  seen  to 
cause  a  noticeable  impact  on  some  high-speed  flows,  such  as 
micronozzle  expansions.7 

The  applicability  of  continuum  numerical  approaches, 
such  as  those  based  on  the  solution  of  Navier-Stokes 
equations,  is  limited  for  flows  where  the  roughness  size  is 
smaller  than  the  gas  mean  free  path.  There  has  been  work 
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recently  on  slip  conditions  which  overcome  this  limitation.4 
In  most  of  the  aforementioned  applications,  the  characteristic 
roughness  size  is  smaller  than  1  /im,  and  is  usually  on  the 
order  of  or  smaller  than  the  gas  mean  free  path.  Being  essen¬ 
tially  a  kinetic  effect,  the  interaction  of  gas  molecules  with 
such  small  surface  roughness  requires  a  kinetic  numerical 
approach  in  order  to  properly  address  the  gas-surface  inter¬ 
action  at  the  level  of  molecular  velocity  distribution  func¬ 
tions.  Kinetic  modeling  of  surface  roughness  effects10  is  hin¬ 
dered  by  a  number  of  problems,  from  the  accurate 
representation  of  complex  shapes  and  roughness  patterns  to 
difficulties  of  experimental  validation. 

The  former  has  been  approached  by  a  number  of  models, 
different  in  their  numerical  complexity  and  physical  accu¬ 
racy.  Among  the  simplified  models  are  the  Cercignani- 
Lampis  (CL)  kernel11'12  and  cone  model.1314  In  the  CL 
model,  rough  surfaces  may  be  simulated  by  specifying  the 
tangential  accommodation  coefficient  a,  larger  than  unity.15 
A  similar  idea  was  applied  by  Ketsdever  et  al.16  where  an 
antispecular  reflection  was  introduced.  In  that  model,  a  fixed 
fraction  of  molecules,  a,  is  reflected  so  that  the  reflected 
normal  and  tangential  momentum  is  the  same  magnitude  as 
the  incoming,  but  with  opposite  directions.  The  models1116 
are  easy  to  implement  and  may  be  used  with  particle  as  well 
as  deterministic  kinetic  solvers,  but  they  do  not  have  a  clear 
correlation  between  the  parameters  of  the  models  and  the 
physical  properties  of  the  rough  surface.  This  drawback  was 
overcome  for  particle  simulations  in  the  cone  models1314 
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where  the  surface  is  either  represented  by  a  number  of  virtual 
cones  or  precomputed  distribution  functions  are  used  at  the 
processing  stage.  These  models  refer  to  the  actual  average 
angle  of  the  rough  surface,  although  it  is  hard  to  use  a  pre¬ 
scribed  angular  distribution,  and  a  special  cut-off  has  to  be 
used  limiting  the  number  of  cones.  Another  approach  is  to 
explicitly  specify  the  rough  surface  as  a  set  of  randomized  or 
regular  structures.  While  this  is  a  reasonable  approach  for 
simple  geometries,  it  may  be  difficult  to  use  for  complex, 
especially  three-dimensional,  surfaces  and  roughness  struc¬ 
ture.  Among  the  more  sophisticated,  although  difficult  to 
implement,  roughness  models  is  the  fractal  model18'14  and  an 
analytic  model.”0 

The  experimental  validation  of  rough  surface  models  is 
hampered  by  requirements  of  high  precision  and  small  scale, 
both  of  which  increase  the  difficulty  of  experimental  mea¬ 
surements,  as  well  as  the  superposition  of  different  physical 
effects.”1  Examples  of  related  work  on  surface  roughness 
are18  where  flow  rates  through  short  channels  were  measured 
and  calculated  with  a  Monte  Carlo  method,  and,”''  where 
triangular  roughness  in  channel  flows  was  studied  numeri¬ 
cally  and  experimentally  for  large  Knudsen  numbers.  Adding 
to  the  difficulty  of  numerical  validation  is  the  complexity  of 
relating  the  measured  macroparameters  to  the  microphenom¬ 
enon  of  interest  to  the  investigation.  Although  changes  in 
macroparameters,  such  as  pressure  distribution  and  mass 
flow,  can  be  affected  by  changes  in  microparameters,  such  as 
accommodation  coefficient,  these  macroparameters  do  not 
necessarily  indicate  specific  reasons  for  the  change  in  accom- 
modation.  Changes  in  angle  of  incidence,  gas  temperature, 
surface  temperature,  surface  material,  and  surface  cleanliness 
can  all  lead  to  changes  in  accommodation.  To  complicate 
matters,  these  parameters  can  be  difficult  to  quantify  experi¬ 
mentally.  In  addition  to  the  change  in  accommodation  due  to 
the  surface,  a  long  microchannel  flow  may  be  affected  by 
variation  of  accommodation  with  local  gas  pressure.  In  most 
cases,  the  probability  of  surface  sticking/accommodation  de¬ 
creases  as  the  availability  of  surface  sites  decrease. 
Lundstrom”4  observed  the  surprising  result  that  Knudsen 
flow  diffusion  increased  with  backing  pressure.  This  was  hy¬ 
pothesized  to  be  due  to  the  variation  in  sticking  with  local 
gas  density  and  the  tendency  of  molecules  to  scatter  more 
specularly  off  an  absorbed  wall  gas  than  off  the  surface. 

This  study  consists  of  a  numerical  investigation  and  ex¬ 
perimental  validation  of  the  effects  of  surface  roughness. 
Specifically,  the  scope  will  be  limited  to  two-dimensional 
helium  and  nitrogen  microchannel  flows  with  well- 
characterized  surface  roughness.  These  flows  will  be  in  the 
low  Reynolds  number  regime,  based  on  entrance  conditions, 
with  surface  rarefaction  from  near  free  molecular  to  near 
continuum,  based  on  surface  roughness.  A  simple  model  of 
the  surface  roughness  that  maintains  the  detailed  balance  and 
is  applicable  to  the  direct  simulation  Monte  Carlo  (DSMC) 
method  has  been  developed,  and  its  connection  to  the 
Cercignani-Lampis  scattering  kernel  is  established.  This  new 
model  is  more  efficient  in  both  setup  cost  and  simulation 
time  than  the  direct  representation  of  the  physical  geometry. 
It  does  not  exhibit  the  shortcomings  seen  in  previous  simple 
surface  roughness  models  while  being  versatile  enough  to 


FIG.  1.  Schematic  of  new  conical  model  reflection  algorithm. 

incorporate  either  an  average  surface  angle  or  a  distribution 
of  surface  angles,  as  well  as  arbitrary  roughness  coverage. 

II.  NUMERICAL  APPROACH 
A.  Conical  model 

Molecular  collisions  with  a  rough  surface  can  be  sepa¬ 
rated  into  two  stages:  (1)  a  molecule  strikes  the  wall  at  a 
particular  point  and  (2)  the  molecule  experiences  one  or 
more  collisions  with  the  surface  roughness  feature  before 
leaving  the  surface.  When  the  characteristic  roughness  size  is 
significantly  less  than  the  mean  free  path  of  the  gas,  no  in- 
termolecular  collisions  occur  during  the  second  stage.  The 
time  taken  for  this  stage  is  much  smaller  than  the  mean  col¬ 
lision  time,  and  the  distance  traveled  is  on  the  order  of  the 
roughness  size.  Therefore,  for  numerical  simulation,  the  time 
and  distance  can  be  ignored  and  it  is  sufficient  to  specify 
only  the  reflected  velocity  of  the  molecule  after  the  collision. 

To  represent  this  two-stage  process  in  an  efficient,  yet 
representative  way,  a  conical  model  is  used  to  select  the  re¬ 
flected  velocity  of  the  collided  molecule.  At  the  point  of 
collision,  the  rough  surface  is  represented  as  a  randomly  po¬ 
sitioned  virtual  conical  hole  with  fixed  opening  angle  /3  and 
height  h.  The  molecule  is  assumed  to  enter  the  cone  through 
the  base  at  point  A,  as  seen  in  Fig.  1.  The  collision  point  for 
the  molecule  and  the  cone  wall  is  then  calculated  at  point  B. 
Here  the  molecule’s  reflected  velocity  is  selected  assuming  a 
diffuse  reflection,  although  a  different  reflection  model  can 
be  used.  The  process  of  calculating  subsequent  collision 
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FIG.  2.  Distribution  function  of  the  angle  between  the  reflected  tangential 
velocity  and  the  flow  direction,  M- 0.3,  /3= 66°. 


points  and  determining  postcollision  velocities  is  continued 
until  the  molecule  leaves  the  cone.  The  point  of  origin  for  the 
reflected  molecule  is  assumed  to  be  the  same  as  the  incoming 
point  of  collision;  all  movement  in  the  virtual  cone  is  only  a 
means  for  determining  the  reflected  velocity  of  the  colliding 
molecule.  Therefore,  the  actual  value  of  h  is  not  important 
(e.g.,  h=  1  m  can  be  used).  In  this  algorithm,  there  are  no 
uncertainties  associated  with  cutoff  for  very  long  traveling 
distances  such  as  those  in  Ref.  13.  Due  to  surface  shape 
simplification,  this  algorithm  requires  only  one  parameter, 
the  average  surface  slope  (represented  as  the  cone  angle  /3), 
making  it  easy  to  implement  in  DSMC. 

To  study  the  influence  of  a  rough  wall,  represented  by 
the  conical  model,  on  flow  properties,  examine  how  the  ve¬ 
locity  distribution  of  incident  molecules  is  transformed  by 
collisions  with  the  wall.  Consider  a  distribution  of  incoming 
molecules  from  an  equilibrium  flow  characterized  by  a  Mach 
number  M  and  directed  along  the  y-axis,  assuming  the  rough 
surface  to  be  in  the  xy  plane.  This  scenario  has  been  simu¬ 
lated  using  the  algorithm  described  above  and  the  distribu¬ 
tion  of  the  angle,  ip,  between  the  reflected  tangential  molecu¬ 
lar  velocity  and  the  flow  direction  is  plotted  in  Fig.  2.  Note 
that  the  maximum  of  the  distribution  corresponds  to  the  di¬ 
rection  opposite  to  the  flow,  which  is  generally  expected  to 
increase  flow  resistance. 

Verification  for  the  conical  surface  model  has  been  per¬ 
formed  for  a  two-dimensional  thermal  bath  to  assure  adher¬ 
ence  to  detailed  balance.  The  thermal  bath  constituted  helium 
gas  initially  heated  to  1000  K  within  a  2D  box  with  a  surface 
temperature  of  300  K.  Gas-surface  collisions  were  simulated 
by  the  conical  roughness  model  with  a  cone  angle  of  (3 
=45°.  The  temperature  relaxation  time  inside  the  test  box  is 
illustrated  in  Fig.  3  for  a  cross  section  along  the  centerline. 
As  particles  collide  with  the  surface,  the  temperature  de¬ 
creases  from  its  initial  value  to  the  equilibrium  value  of 
300  K.  By  coming  to  a  stable  equilibrium,  the  simulation 
proves  adherence  to  the  detailed  principle. 


FIG.  3.  Thermal  relaxation  of  helium  in  a  300  K  box.  Time  is  in  mean 
collision  times,  rx,  at  1000  K. 


B.  Cercignani-Lampis  scattering  kernel 

The  above  stated  cone  model  technique  is  valuable  for 
particle  simulations,  but  does  not  have  a  direct  applicability 
in  deterministic  kinetic  solvers  since  it  is  based  on  individual 
particle  velocities  and  not  on  distribution  functions.  It  is  rea¬ 
sonable  to  establish  a  connection  between  the  conical  model 
and  a  kernel  that  is  directly  applicable  to  a  kinetic  solver. 
One  such  kernel,  the  Cercignani-Lampis  (CL)  (Refs.  11  and 
25)  scattering  kernel,  can  be  used  in  kinetic  solvers  but  lacks 
correlation  between  the  physical  surface  and  the  main  input 
parameter  for  the  kernel.  The  conical  surface  model  can  be 
used  to  relate  this  parameter,  the  tangential  momentum  ac¬ 
commodation  coefficient  a,,  to  an  average  surface  roughness 
angle. 

To  facilitate  this  relation  it  is  necessary  to  create  an  ana¬ 
lytic  expression  for  the  distribution  of  <p  as  a  function  of  a,. 
Assuming  the  temperature  of  the  gas  and  the  wall  are  equal 
and  velocities  are  normalized  by  ym/kT ,  the  distribution 
functions  of  x  and  y  velocity  components  of  incident  mol¬ 
ecules  are 


jtw = 


exp(-  (v  -M\  y)2/2), 


(1) 


-f=  exp(-  tr/2), 
\I2tt 


(2) 


respectively.  After  the  CL  transformation,  the  x  component 
of  the  reflected  velocity  is  equal  to 


ref / 


1  at  inc  CL 
,  Vx  +V x  ’ 

V 1  -  a, 


(3) 


where  the  distribution  function  of  i\L  is  given  by 


1 


V27t(1  -  |1  -  a, |) 


exp  - 


2tt(1  -|1  -  a, \) 


(4) 


Here,  a,  is  the  CL  parameter  that  corresponds  to  the  tangen¬ 
tial  accommodation  coefficient,  0Sfff<2.  a,  =  I  corresponds 
to  zero  average  tangential  momentum  relative  to  the  surface; 
the  larger  the  deviation  between  a,  and  unity,  the  larger  the 
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FIG.  4.  Tangential  momentum  accommodation  coefficient,  at,  as  a  function 
of  cone  opening  angle  (3,  for  M=03. 


average  momentum  of  the  reflected  molecules.  Note  that  v'f1 
is  the  sum  of  the  two  independent  normally  distributed  com¬ 
ponents,  therefore,  its  distribution  function  can  be  written  as 

=  ).  (5) 

The  distribution  function  of  the  y  velocity  component  does 
not  change  during  the  CL  transformation.  Note  that  tan(  </>) 
=  |uf'|/^ef/.  so  f°r  0<(/><  7r/2  the  distribution  function  of 
<b  can  be  obtained  with 

roo 

/(tan  4>)  =  I  u/f{u^ff(ux  tan  (f>)dux,  (6) 

Jo 

and  for  tt/2  <cf><TT 

/(tan  (/>)  =  -[  ujrf‘(.ux)fycn(ux  tan  (f>)du  x .  (7) 

J  —oo 

Finally,  the  reflected  velocity  distribution  function  can  be 
written 

M  4>)  =  2/(tan  4>)  ( 1  +  tan2  cf>)  (8) 

or 

-f2  JL 

e  t  cos  <p  . 

/,/,(/')  = - +  — T=~  exp(-  t  sin-  </)(l  +  erf(/  cos  </)), 

77  Vtt 

(9) 

where  t=My  y{  1  -  a,)  /  V 1 1  -  a, \ . 

Parameter  a,  of  the  CL  transformation  can  be  found  by 
least  square  fitting  of  </>  distributions  obtained  using  the  cone 
reflection  algorithm  described  in  Sec.  II  A  as  seen  in  Fig.  2. 
Figure  4  shows  the  value  of  a„  found  by  a  least  square  fit, 
for  M=0.3  and  different  cone  angles  (i.  As  expected,  a ,  is 
maximum  at  some  intermediate  value  of  the  cone  angle. 
Since  a,  is  larger  than  1,  the  average  tangential  momentum 


FIG.  5.  Tangential  momentum  accommodation  coefficient,  at  as  a  function 
of  Mach  number  for  two  cone  opening  angles. 


of  the  reflected  molecules  points  in  the  opposite  direction 
with  respect  to  the  average  tangential  momentum  of  the  in¬ 
coming  molecules. 

The  value  of  a,  only  weakly  depends  on  Mach  number, 
as  seen  in  Fig.  5,  where  a,  is  shown  as  a  function  of  Mach 
number  for  two  different  opening  angles,  /3=120°  and  /3 
=  66°.  This  fact  facilitates  the  use  of  the  CL  model  in  simu¬ 
lations,  since  the  value  of  a,  can  be  selected  universally  de¬ 
pending  on  the  degree  of  roughness  of  the  surface,  and  not 
on  the  flow  properties. 

C.  Simulation  conditions 

In  order  to  validate  the  conical  roughness  model  and  its 
link  to  the  CL  scattering  kernel,  both  were  individually  in¬ 
corporated  into  a  DSMC  solver  and  simulations  were  run  for 
comparison  against  experiment.  To  maximize  the  measurable 
effect  of  the  different  surface  models  on  the  flow,  the  simu¬ 
lations  were  conducted  on  a  long,  narrow  microchannel  with 
an  LIH=  100  and  an  assumed  infinite  width.  The  2D  capa¬ 
bility  of  SMILE  (Ref.  26)  was  used  in  all  DSMC  computa¬ 
tions.  This  system  utilizes  the  majorant  frequency  scheme  to 
calculate  intermolecular  interactions,  while  using  the  vari¬ 
able  hard  sphere  model  for  intermolecular  potentials.  Four 
surfaces  were  used:  (i)  a  fully  diffuse  smooth  surface  with 
complete  energy  and  momentum  accommodation,  (ii)  a 
smooth  surface  utilizing  the  conical  surface  roughness  model 
with  an  opening  angle  of  66°  which  corresponded  to  an  ear¬ 
lier  experimental  setup  concept,  (iii)  CL  model  with  a  tan¬ 
gential  momentum  accommodation  coefficient  of  1.045  at¬ 
tained  using  the  conical  model  with  a  66°  cone  opening 
angle,  and  (iv)  a  diffuse  surface  consisting  of  about  900  tri¬ 
angles  representative  of  an  actual  experimental  shape. 

The  channel  length  and  half  height,  from  plane  of  sym¬ 
metry  to  the  tops  of  the  surface  features,  were  1.5  cm  and 
75  /i,m,  respectively.  Helium  was  used  as  the  test  gas  with 
stagnation  pressure  from  200  Pa  to  8000  Pa  for  all  four  sur¬ 
faces.  The  Knudsen  number  is  defined  as 
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FIG.  6.  Schematic  of  the  computational  domain. 


(b)  Flow 


X  1 _ kT 

L  V 2mrcl2L  yflPircfL 


(10) 


where  P  is  the  stagnation  pressure,  k  is  Boltzmann’s  con¬ 
stant,  T  is  the  stagnation  temperature  (300  K),  d  is  the  mo¬ 
lecular  diameter,  and  L  is  the  feature  size  (about  20  /tm). 
The  Knudsen  number  for  these  tests  ranged  from  0.1  to  5. 
Nitrogen  was  used  as  the  test  gas  for  only  surfaces  (i)  and 
(ii),  with  stagnation  pressures  from  200  Pa  to  6000  Pa. 
Again,  based  on  the  stagnation  conditions  and  the  surface 
roughness  size,  the  Knudsen  number  ranged  from  0.05  to  2. 
A  convergence  study  was  performed  to  assess  the  impact  of 
the  inflow  boundary  conditions  on  the  results  of  the  compu¬ 
tations.  These  computations  were  performed  using  helium  at 
6000  Pa  and  two  inflow  boundary  locations:  the  first  extend¬ 
ing  1700  /tm  from  the  channel  entrance  and  1700  /xm  from 
the  symmetry  plane,  the  other  850  and  850  /tm,  respectively. 
A  Maxwellian  distribution  function  with  zero  flow  velocity 
was  assumed  at  these  boundaries.  No  impact  from  the 
boundary  location  was  found.  The  larger  domain  was  used  in 
all  subsequent  computations,  with  the  number  of  molecules 
ranging  from  about  1.6  million  for  the  lower  pressures  to 
about  10  million  for  the  larger  ones.  A  schematic  for  the 
computational  domain  can  be  seen  in  Fig.  6,  where  the  out¬ 
flow  boundary  condition  is  modeled  as  a  vacuum. 


III.  EXPERIMENTAL  APPROACH 

The  experiments  were  conducted  using  atomically 
smooth,  triangular  patterned,  and  square  patterned  silicon 
textures.  Standard  microelectromechanical  systems  (MEMS) 
processes  were  performed  to  create  1  cm  by  1.5  cm  textured 
areas  centered  on  2  cm  by  1.5  cm  chips.  The  feature  width 
for  both  square  and  triangular  textures  is  approximately 
15  /jl m.  The  triangles  are  approximately  11  /im  deep  with  an 
etched  angle  of  54.7°  from  surface  parallel,  which  would 
give  an  opening  angle  of  approximately  70°,  very  close  to 
the  numerical  simulation.  The  squares  are  approximately 
30  /im  deep,  and  both  geometries  can  be  seen  in  Fig.  7. 
Teflon  shims  were  placed  between  the  chips  to  create  a  chan¬ 
nel  1  cm  wide,  150  /jl m  high,  and  1.5  cm  long  such  that  gas 
flow  along  the  length  of  the  channel  was  perpendicular  to  the 
textures  etched  in  the  silicon  chip.  Along  this  length  there 
were  approximately  900  features.  This  channel  was  secured 
in  an  aluminum  holder  and  dimensions  of  the  channel  were 
measured  with  a  scanning  electron  microscope  before  and 
after  testing  as  seen  in  Fig.  8.  The  ratio  of  the  measured  open 
area  of  the  channel  to  the  open  area  of  the  numerical  simu¬ 
lation  was  used  to  linearly  correct  the  experimental  mass 


FIG.  7.  Average  experimental  surface  roughness  feature  dimensions  for  (a) 
the  triangular  surface  and  (b)  the  square  surface. 


flow  data.  Typical  deviation  for  the  channel  height  was  3%- 
5%.  Deviation  for  the  width  was  less  than  1%. 

The  aluminum-Teflon-silicon  assembly  was  placed  in  the 
wall  of  a  3550  cm3  plenum,  which  acted  as  a  stagnation 


FIG.  8.  SEM  pictures  of  (a)  the  end  of  the  triangular  textures  (from  above), 
and  (b)  the  view  along  cleaved  square  textures  (from  side). 
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FIG.  9.  Normalized  pressure  and  Mach  number  as  a  function  of  distance 
downstream  along  the  centerline  for  two  stagnation  conditions. 

chamber.  The  plenum  was  machined  with  two  gas  inlet  ports 
and  three  ports  for  pressure  measurements.  The  setup  was 
operated  in  the  CHAFF  IV  vacuum  chamber  at  the  Univer¬ 
sity  of  Southern  California,  which  is  capable  of  an  ultimate 
pressure  of  10-6  Torr  with  working  pressures  no  higher  than 
10-4  Torr.  This  allowed  plenum  pressures  from  approxi¬ 
mately  1  to  100  Torr  with  ratios  of  stagnation  to  ambient 
pressures  greater  than  10s  maintained  throughout.  Helium  or 
nitrogen  test  gases  were  used.  Upstream  gas  flow  was  mea¬ 
sured  with  a  mass  flow  meter  and  regulated  through  a  needle 
valve.  The  pressure  in  the  plenum  was  monitored  with  sev¬ 
eral  differential  baratrons.  This  array  of  instrumentation  al¬ 
lowed  measurements  of  pressure  and  mass  flow  with  instru¬ 
mental  precisions  below  1%. 

IV.  RESULTS  AND  DISCUSSION 

The  first  result  to  consider  is  that  of  the  numerical  simu¬ 
lation  using  the  smooth,  fully  diffuse  surface  with  a  helium 
test  gas.  Figure  9  shows  the  local  pressure,  normalized  by  the 
stagnation  value,  and  Mach  number  along  the  centerline  of 
the  channel  for  plenum  pressures  of  200  Pa  and  2000  Pa. 
The  Knudsen  number,  based  on  the  channel  height  of 
150  /J.m  instead  of  the  feature  size,  for  the  200  Pa  case  in¬ 
creases  from  below  1  to  about  100.  As  seen  in  Fig.  9,  the 
pressure  in  the  free  molecular  flow  drops  linearly  from  P0 
=  200  Pa  to  0.0 1  P{]  through  the  channel.  The  Mach  number 
stays  below  0.05  for  the  majority  of  the  channel,  increasing 
to  1  at  the  exit;  the  Mach  number  for  the  P0  =  2000  Pa  case  is 
higher,  as  expected.  This  allows  for  the  use  of  a  constant  a, 
value  in  the  CL  model  along  the  length  of  the  channel. 

Figure  10  shows  the  ratio  of  local  pressure  for  the  coni¬ 
cal  surface  simulation  divided  by  the  local  pressure  for  the 
smooth  surface  as  a  function  of  longitudinal  distance  down 
the  channel.  As  expected,  the  pressure  along  the  channel  in¬ 
creases  for  the  rough  surface  compared  to  the  smooth,  since 
the  rough  surface  impedes  the  flow.  This  effect  decreases  as 
the  flow  becomes  free  molecular,  at  which  point  the  pressure 
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FIG.  10.  Pressure  ratio  as  a  function  of  distance  downstream  along  the 
centerline,  conical  rough  surface  to  the  smooth  surface. 

ratio  down  the  channel  asymptotes  to  1  far  from  the  exit.  The 
ratio  drops  below  1  at  the  point  where  the  local  mean  free 
path  exceeds  the  channel  height  allowing  molecules  to  “see” 
the  exit  without  intermolecular  collisions.  Downstream  from 
this  point,  any  molecules  traveling  downstream  are  more 
likely  to  leave  the  channel  after  colliding  with  the  rough 
surface  than  they  would  be  after  colliding  with  a  smooth 
surface.  Upstream  from  this  point,  molecular  collisions  are 
likely  to  prevent  this.  Note  that  a  relatively  small  difference 
between  pressures  for  the  rough  and  smooth  surfaces  is  ob¬ 
served  for  much  of  the  channel  length,  consistent  with  the 
measurements  of  Turner  et  al.6 

Calculated  helium  mass  flow  as  a  function  of  stagnation 
pressure  is  listed  in  Table  I  for  all  four  surfaces.  A  reduction 
in  mass  flow  of  25%  at  the  lower  pressures  and  6%  at  higher 
pressures  is  seen  between  the  conical  surface  model  and  the 
smooth,  fully  diffuse  case.  The  CL  model  with  an  a, 
=  1.045  matches  the  conical  model  to  7%  at  the  lower  pres¬ 
sures  and  less  than  1%  at  the  higher  pressures.  The  conical 
surface  matches  the  triangular  geometry  pattern  within  3%  at 
the  lower  pressures  and  less  than  1%  at  the  higher  pressures. 
The  impact  of  the  surface  roughness  is  larger  for  helium  than 
for  nitrogen  at  any  given  pressure;  however,  when  plotted 
versus  Knudsen  number,  the  impact  is  similar.  Assumptions 
made  for  the  cone  surface  model  require  a  noncontinuum 


TABLE  I.  Numerical  results  for  mass  flow  as  a  function  of  stagnation  pres¬ 
sure  for  various  surface  models  using  a  helium  test  gas. 


Pressure 

Pa 

Surface  (i) 
kg/s 

Surface  (ii) 
kg/s 

Surface  (iii) 
kg/s 

Surface  (iv) 
kg/s 

200 

4.132E-09 

3.188E-09 

3.440E-09 

3.093E-09 

650 

1.381E-08 

1.125E-08 

1.162E-08 

1.105E-08 

2000 

5.680E-08 

4.939E-08 

4.997E-08 

4.898E-08 

4000 

1.623E-07 

1.484E-07 

1.502E-07 

1.483E-07 

6000 

3.094E-07 

2.899E-07 

2.893E-07 

2.896E-07 

8000 

5.098E-07 

4.810E-07 

4.830E-07 

4.810E-07 
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FIG.  11.  Comparison  of  the  experimental  smooth  texture  geometry  and 
numerical  smooth  diffuse  roughness  model  for  (a)  helium  gas,  (b)  nitrogen 
gas. 
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FIG.  12.  Comparison  of  the  experimental  triangular  texture  geometry  and 
numerical  cone  roughness  model  for  (a)  helium  gas,  (b)  nitrogen  gas. 


Knudsen  number  based  on  the  feature  size  of  the  roughness. 
Although  the  cone  model  agrees  favorably  with  the  direct 
representation  over  the  full  range  tested,  it  is  not  a  given 
conclusion  that  the  cone  model  is  applicable  in  low  Knudsen 
number  cases. 

There  is  good  agreement  between  the  experimental 
smooth  silicon  surface  and  the  smooth  fully  diffuse  simula¬ 
tion  as  seen  in  Fig.  11.  The  experimental  uncertainty,  based 
on  standard  deviation,  varies  from  about  10%  at  P0 
<200  Pa  to  less  than  1%  at  PQ>  10000  Pa.  There  is  a  2% 
underprediction  by  the  DSMC  simulations  along  the  range  of 
pressures  tested.  One  of  three  effects  can  cause  this  devia¬ 
tion.  The  first  is  a  finite  specularity  that  may  have  existed  in 
the  experiment  and  was  unaccounted  for  in  the  simulations. 
The  second  is  an  uncertainty  of  approximately  2%  in  the  area 
of  the  channel  based  on  the  SEM  pictures.  The  third  is  the 
previously  mentioned  influence  of  surface  absorption  and  the 
decrease  of  surface  sites  at  higher  pressures,  also  creating  a 
higher  physical  specularity  then  accounted  for  in  the  simula¬ 
tion.  However,  relative  influence  of  these  effects  are  un¬ 
known.  It  is  important  to  note  that  the  pressure  dependence 
of  the  underprediction  is  hard  to  determine  due  to  the  larger 
experimental  uncertainty  at  lower  pressures. 

There  is  also  good  agreement  between  the  experimental 
triangular  surface  and  the  cone  model  simulations  as  seen  in 
Fig.  12.  As  with  the  smooth  surface,  the  experimental  uncer¬ 
tainty,  based  on  standard  deviation,  varies  from  about  10%  at 
P0<200  Pa  to  less  than  1%  at  P0>  10000  Pa.  At  the  higher 
pressures,  there  is  less  than  2%  difference  between  the  ex¬ 


periment  and  simulation;  however,  at  the  lower  pressures, 
there  is  a  5%  overprediction  by  the  DSMC.  Again  both  the 
uncertainty  in  the  experimental  channel  area  and  the  larger 
statistical  scatter  at  lower  pressures  could  contribute  to  this 
difference. 

The  experimental  data  for  the  square  surface  features  can 
be  seen  in  Fig.  13  compared  with  the  other  two  experimental 
surfaces.  Molecular  interaction  with  the  vertical  walls  of  the 
square  surfaces  could  be  approximated  as  cones  with  open¬ 
ing  angles  (3=0 °.  As  expected,  the  square  data  falls  between 
the  smooth  (3=  1 80°  and  triangular  (3=10°  experimental 
points.  Although  the  square  surface  data  trends  towards  the 
triangular  data  at  lower  pressures,  it  moves  towards  the 
smooth  data  at  higher  pressures.  This  could  be  attributed  to 
the  influence  of  the  exit  plane  on  the  flow  and  the  reduction 
in  Knudsen  number  as  the  pressure  ratio  increases.  At  higher 
pressures  and  densities,  the  square  texture  acts  more  like  a 
fully  diffuse  surface  than  a  series  of  roughness  features,  until 
a  few  mean  free  paths  before  the  exit  plane,  ultimately  as- 
ymptoting  at  the  smooth  surface  value. 

V.  CONCLUSION 

A  conical  surface  roughness  model  has  been  developed, 
numerically  verified,  and  experimentally  validated.  For  a 
particle  simulation,  the  conical  model  reduces  both  setup  and 
simulation  time  over  direct  geometric  representation  without 
sacrificing  accuracy.  For  deterministic  solvers,  this  conical 
model  gives  a  crucial  link  between  the  Cercignani-Lampis 
(CL)  scattering  kernel,  which  requires  a  defined  tangential 


Downloaded  10  Oct  2007  to  198.14.54.253.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pof.aip.org/pof/copyright.jsp 


106101-8  Lilly  et  al. 


0  2000  4000  6000  8000  10000  12000  14000  16000 

(b)  Pressure  (Pa) 


FIG.  13.  Experimental  results  for  the  textured  geometries  for  helium  gas. 
(a)  is  a  zoom  of  the  lower  portion  of  (b). 

momentum  accommodation  coefficient,  a„  and  a  physical 
surface  characteristic,  the  average  surface  roughness  angle. 
The  conical  model  is  general  enough  to  be  used  for  surfaces 
with  both  prescribed  roughness  features  and  nonuniform 
roughness  features.  In  the  latter  case,  a  mix  of  conical  open¬ 
ing  angles  and  traditional  diffuse  interactions  can  be  used 
based  on  the  surface  distribution. 

Numerical  simulations  using  the  conical  model,  the  CL 
scattering  kernel,  and  a  direct  geometric  representation  of  a 
regular  triangular  surface  were  compared  against  a  smooth, 
fully  diffuse  model.  All  four  cases  were  simulated  with  the 
DSMC  method  as  two-dimensional  flows  of  helium  or  nitro¬ 
gen  through  a  150  /jl m  tall,  1.5  cm  long  microchannel  with 
Reynolds  numbers  from  0.01  to  10  based  on  inlet  conditions. 
The  direct  geometric  representation,  the  most  costly  but  most 
accurate  compared  to  experiment,  showed  a  decrease  in  mass 
flow  between  25%  for  the  lower  stagnation  pressures  simu¬ 
lated  and  6%  for  the  higher  stagnation  pressures  simulated. 
The  conical  model  matched  this  result  to  within  3%  at  the 
lower  stagnation  pressures  and  less  than  1%  for  the  higher 
stagnation  pressures  while  the  CL  model  matched  this  result 
within  7%  and  less  than  1%,  respectively. 

Experiments  were  conducted  on  transitional  flows, 
where  the  mean  free  path  is  on  the  order  of  or  larger  than  the 
roughness  size,  through  a  1  cm  wide  channel  with  similar 
height  and  length  to  the  simulations  to  approximate  2D  con¬ 
ditions.  The  channel  walls  were  made  of  silicon  with:  (i) 
polished  smooth  surfaces;  (ii)  regular  triangular  roughness; 
and  (iii)  regular  square  roughness  with  characteristic  rough¬ 
ness  scales  of  <  1  /jl m,  1 1  /im,  and  29  /xrn,  respectively.  Ex- 
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perimental  results  validate  the  numerical  simulations,  which 
match  to  within  2%  throughout  the  range  tested  for  the 
smooth  geometry  and  between  2%  and  5%  for  the  high  and 
low  stagnation  pressures,  respectively,  for  the  triangular 
roughness  geometry.  The  experimental  mass  flow  for  the 
square  roughness  geometry  falls  between  the  smooth  and 
triangular  cases  over  the  entire  test  range  as  expected  based 
on  yS. 
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